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MissourRI AGRICULTURAL EXPERIMENT STATION 


Communicated April 13, 1944 


No dominant mutations have been found in x-ray progenies of barley and 
maize, in the course of experiments in which many hundreds of recessive 
mutations have been detected. Many dominant alterations appear in the 
first generation progeny produced by the use of x-rayed pollen, but these 
alterations are lethal to the gametophyte and are therefore not reproduced 
in later generations. The recessive mutations detected are at least par- 
tially haplo-viable, since they must be transmitted through both male and 
female gametophytes before their effects may be observed in homozygous 
individuals. Since all of the induced alterations which are viable in haplo- 
phase are wholly recessive, these experiments provide no evidence of the 
production by x-rays of any mutant gene of positive action. The induced 
mutations observed may be merely gene losses tolerated by the gameto- 
phyte.?? 

The evidence against the occurrence of x-ray-induced dominant muta- 
tion is, however, inconclusive, for the following reasons: 

1. The number of genes capable of showing the effects of dominant 
mutation may be much smaller than the number capable of showing re- 
cessive mutation, since many genes may be fixed by natural selection at a 
level maximal for phenotypic effect. The preponderance of recessive mu- 
tations may be only a reflection of the preponderance of such alleles in the 
genotype of the irradiated individual. The possibility of inducing domi- 
nant mutation therefore may be tested critically only by the effects upon 
known recessives, that is, genes with known dominant alleles. 

2. Among known recessives many may be themselves deficiencies and 
therefore incapable of mutation. Critical evidence of failure to mutate to 
a dominant allele therefore may be obtained only from recessive genes 
which have previously been known to mutate to dominant alleles. 

3. The only recessive alleles which meet this requirement are the so- 
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called variegation genes, which may be regarded as unstable recessives 
mutating frequently to a dominant allele. Each of these is characterized 
by a high spontaneous frequency of dominant mutation, so high that the 
effect of x-ray treatment in inducing additional dominant mutations prob- 
ably would be inappreciable. Demerec’s experiment‘ with the unstable 
miniature-3 of Drosophila virilts is an example of this difficulty. 

It is possible to avoid these difficulties in the case of one gene. The re- 
cessive gene a; in maize has several known dominant alleles. Rhoades*® 
has shown that in the presence of another gene, Dt, mutations of a to various 
dominant alleles occur. With increasing dosage of Dt the frequency of 
dominant mutation of a is sharply increased, and in endosperms triplex for 
Dt and for a (a a a Dt Dt Dt) more than 100 mutations per seed are usually 
found. In the absence of Dit no dominant mutations of a have been re- 
ported, though a aa dt dt dt endosperms have been closely examined in 
large numbers. (Dr. Rhoades informs me that he has seen one instance 
of an apparent “‘dot’”’ of anthocyanin colored tissue in the aleurone cells of 
an aaa dt dt dt seed.) 

The gene a, then, is capable of mutation to a recognizable dominant allele, 
and can be tested for response to x-ray treatment in genotypes in which 
spontaneous dominant mutations of the gene rarely if ever occur. Since 
the effects of mutation are clearly recognizable in minute sectors of tissue, 
the treatment may be applied at a fairly advanced stage in endosperm de- 
velopment, so that many hundreds of cells may be tested for mutation by 
examination of a single endosperm. It is possible therefore to test for the 
occurrence of this mutation in practically unlimited populations. 

The experiment was conducted as follows: The seed to be irradiated was 
produced by the cross aa X A a, both parent stocks being homozygous 
for dt dt and for the complementary factors required for aleurone color de- 
velopment. The endosperms of half of the seeds produced are A a a 
(colored aleurone). These serve to indicate the size of the sectors resulting 
from genetic alterations induced by irradiation at the stage of development 
chosen, since induced deficiencies or recessive mutations of A in these en- 
dosperms result in sectors deficient in aleurone color. In the colorless 
seeds (a a a) induced dominant mutation of any one of the a genes would 
result in a corresponding sector of colored aleurone. The colored seeds 
thus provide a basis for estimating the number of opportunities for the 
occurrence of detectable mutation in the colorless seeds, and for compari- 
son of the relative frequency of induced dominant mutation and induced 
loss of A due to recessive mutation or deficiency. 

X-ray treatment was applied to the ears 73-81 hours after pollination, 
following preliminary experiments which indicated that treatment at this 
stage results in sectors of suitable size. Treatments were made on field- 
grown plants by means of a mobile x-ray unit with self-rectifying tube, 
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operated at 88 K.V.P., 4.5 ma. tube current, and 15 inches target distance, 
with no filtration other than that of the overlying tissue. Under these 
conditions the time of treatment is 2.48 minutes per 100 r. The doses 
applied, uncorrected for filtration, were 800 r and 1600 r. 

The lower dose was applied to 21 ears, which yielded a total of 8003 seeds. 
The effect of the treatment was clearly greater on the exposed side of the 
ear than on the opposite side, as indicated by the frequency of deficiency 
sectors and also by the occurrence of “‘scarred’’ endosperms. There was 
also appreciable variation in dosage along the long axis. 

The average size and frequency of sectors showing ‘‘A-loss’’ (that is, 
loss of the A phenotype), was estimated by detailed examination of a small 
sample of seeds. From each of six representative ears given the lower dose 
of radiation, three seeds were taken from the treated side and three from 
the opposite side, with approximately even distribution along the long axis 
of the ear. Each of the 36 seeds was closely examined under magnification 
by the use of a dissecting microscope. This permits determination of the 
size of sectors by direct count of aleurone cells affected. In small sectors 
due to loss of A the aleurone cells are not wholly colorless, but the contrast 
in color is sharp enough to make the outline of the sector quite definite. In 
large sectors, produced by irradiation of the ear about 30 hours after poll- 
ination, the aleurone cells of sectors resulting from loss of A are colorless 
except for a rim of faintly colored cells about the margin of the sector. 
There are in addition occasional sectors of dilute aleurone color, probably 
resulting from similar changes affecting an unknown gene or genes. These 
are readily distinguishable from sectors involving A in the early-treated 
material, but are less distinct when the sectors are so small that marginal 
color may cover the entire sector. In order to minimize the error from 
this source, all sectors which seemed possibly due to loss of genes other 
than A were omitted. These made up nearly one-third of the aleurone 
sectors observed on the treated seeds. The ratio of dilute sectors to color- 
less sectors on early-treated seeds of the same cross is considerably less than 
1:2, and it is therefore probable that the actual frequency of A-loss in the 
experimental material is somewhat higher than the data below would indi- 
cate. The observations are summarized in table 1. 

The 18 seeds representing the exposed side of the treated ears (seeds Nos. 
1 to 3 on each ear) yielded 618 sectors, or an average of 34 per seed; while 
those representing the opposite side yielded 248, or an average of 14 per 
seed. The sectors varied rather widely in size, the average number of 
aleurone cells included being 112. The frequency of A-loss sectors on un- 
treated ears, determined by similar examination of 10 representative seeds 
from each of the 19 control ears, was 0.12 per seed. The 23 sectors found 
ranged in size from 9 to 750 cells. 

The approximate number of aleurone cells in the entire endosperm was 
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determined similarly by counts of aleurone cells in representative areas. 
The number of aleurone cells in seeds of average size was found to be about 
160,000. It may be roughly estimated, therefore, that at the time of 
treatment the number of cells per seed whose genetic alteration could be 
detected by sectors in the mature aleurone is about 1400 (160,000/112). 

The colorless seeds of the same treated ears permit the detection of muta- 
tions of a to any one of the various A alleles producing aleurone color. 
Each endosperm cell contains three a’s, each of which is capable, under the 
influence of Di, of mutating to a colored aleurone allele. If the phenotypic 
effects of mutation and of deficiency are shown by the total cell progeny of 
the irradiated cell, or by the same fraction of this progeny, the result of 
mutation should be a colored sector equal in area to the color-deficient sec- 
tors produced by deficiency. 

If either alteration requires for its realization a period of one or more cell 
generations in excess of the other, the resulting sectors may be of unequal 
size. It is not to be expected, however, that delay in realization will regu- 
larly reduce the size of sector by one-half for each cell generation. The 
sectors observed include only the cells of the aleurone layer which are in- 
cluded in the block of tissue affected; their reduction in size with advancing 
endosperm development depends upon the pattern of development of the 
tissues concerned. Actually, the average size of deficiency sector observed 
in ears treated at periods after pollination ranging from 73 to 81 hours was 
not changed enough to show any consistent relation among the 21 ears 
treated. Ears treated 96 hours after pollination yielded smaller but 
readily recognizable sectors. It is evident, therefore, that even if there 
were a delay in realization of gene mutation for several cell generations after 
irradiation, the resultant sectors would be detectable. Colored sectors as 
small as 10 cells in size are clearly recognizable on close examination of the 
seeds without magnification, and a single colored aleurone cell on a color- 
less seed is readily identified under magnification. 

All of the colorless seeds were closely examined without magnification 
for sectors of colored aleurone. Occasional discolored areas of pericarp or 
endosperm were found, but under magnification these were readily dis- 
tinguishable from colored aleurone sectors. No sectors of colored aleurone 
were found, among the 3832 seeds examined. On the basis of 1400 cells 
per seed, each containing 3 a genes capable of dominant mutation, this 
represents approximately 16 million chances for detectable dominant mu- 
tation (3832 X 1400 X 3). The treatment to which these cells were ex- 
posed produced in the adjoining colored seeds approximately 100,000 losses 
of A. Since the number of A’s tested for frequency of A-loss was only one- 
third the number of a’s tested for frequency of dominant mutation, the ex- 
periment shows failure to induce dominant mutation in a trial on a scale 
sufficient to produce about 300,000 A -losses. 
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A similar examination was made of the ears irradiated at 1600 r. At this 
dose the number of A-losses on seeds on the exposed side of the ear is too 
great to permit an accurate count. On the opposite side of the ear the fre- 
quency is very high also, but samples representative of the central portion 
of this region (amounting to approximately one-sixth of the ear) were suit- 
able for examination. In 8 of the ears treated a sample of 6 seeds per ear 
was examined, the sample in each case being taken from the two rows least 
exposed to the radiation. The number of A-losses from these samples 
averaged 20.3 per seed. This is about three times the frequency of A- 


TABLE 1 
FREQUENCY AND SIZE OF SECTORS DuE To A-Loss FOLLOWING X-RAY TREATMENT 
Ear No. 870 No. 876 No. 880 No. 881 No. 882 No. 884 Total 


Interval, pollina- 
tion to irradia- 


tion (hrs.) 75 75.5 80.5 81 81.5 75 
Number of sectors 
Seed No. 1 22 32 26 36 49 40 205 
No. 2 30 31 38 33 34 29 195 
No. 3 49 43 34 23 37 32 218 
No. 4 3 9 13 14 6 15 60 
No. 5 20 12 17 12 14 9 84 
No. 6 11 16 13 18 24 22 104 
Number of aleu- 
rone cells per 
sector (fre- 
quency dis- 
tribution) 
1-10 12 12 16 9 23 18 90 
11-50 38 43 48 42 47 41 259 
51-200 71 64 46 63 82 69 395 
201-500 14 23 27 22 11 19 116 
500+ 0 1 4 0 1 0 6 


losses in seeds similarly located on the ears treated at 800 r, and inspection 
of seeds in other regions of the ears indicates that this is a fair estimate of 
the relative frequency of A-losses from the two doses used. The average 
size of sectors was 113 cells. 

The 15 ears treated at 1600 r produced a total of 2487 colorless seeds. 
These were closely examined for sectors of colored aleurone, and in addition 
a sample of 50 seeds from the most heavily treated region of each was 
minutely examined under magnification. No sectors of colored aleurone 
were found, with the exception of a single doubtful spot possibly including 
3 aleurone cells, which could not be positively identified. The frequency of 
A-losses, producible by the same dose in an equal population of A-genes, 
estimated as before, is about 600,000. 
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The susceptibility to dominant mutation of the a genes used in this ex- 
periment was tested by pollinating each of the parental stocks by a homo- 
zygous stock designated a-dl Dt. The gene a-dl is one of a number of ‘‘dot- 
less’’ a alleles arising by mutation of standard a under the influence of Dt. 
It is indistinguishable from a in phenotypic effect, but differs in showing 
almost no dots due to dominant mutation in the presence of Dt. In endo- 
sperm triplex for a-dl and Dt an occasional small dot may be found, as in 
the case of the similar mutant a* described by Rhoades, but the rate is 
negligible in comparison to that of standard a. In the endosperms of 
a a a-dl dt dt Dt, produced by the crosses under discussion, the dots ob- 
served may be ascribed almost wholly to mutation of one of the two a’s 
under the influence of the single dose of Dt. The mutation frequencies for 
x-ray action and Di action are, of course, not directly comparable, since it 
is not possible to apply the two stimuli for comparable periods or to com- 
parable populations of cells. 

The number of dots per seed in endosperms of a a a-dl dt dt Dt produced 
by crosses on the female parental stock was 2.95, on the colorless seeds of the 
male parental stock 2.1. In each case the dots represent the number of 
mutations produced by two a genes in the presence of a single Dit gene. In 
the a aa seeds irradiated, therefore, the number of mutations to A which 
would be expected in the presence of one dose of Dt is about 4 per seed (2.95 
+ 2.1/2). Rhoades‘ has shown that the frequency of mutation of a to A 
is disproportionately increased by increased dosage of Dt, and that about 
17 times as many mutations are produced in Dt Dt Dt seeds as in Dt dt dt 
seeds. The population of 6319 seeds of a a a, which failed to yield muta- 
tions to A under x-ray treatment, was therefore capable of yielding about 
400,000 mutations to A under the influence of homozygous Dt (6319 X 4.0 
Xx 17). 

Summary.—X-ray treatment failed to induce mutation of a to A (or to 
any other colored aleurone allele), in populations which, with the x-ray 
doses applied, were capable of yielding about 900,000 losses of A by de- 
ficiency or by mutation to a colorless allele. 

The a genes used in the experiment, when combined with Dt, yielded 
numerous mutations to A or other colored aleurone alleles. The population 
irradiated was large enough to have yielded, under the influence of homozy- 
gous Dt, about 400,000 of such mutations. 

1 Coéperative investigations of the Division of Cereal Crops and Diseases, Bureau of 
Plant Industry, Soils, and Agricultural Engineering, Agricultural Research Administra- 
tion, U. S. Department of Agriculture, and Missouri Agricultural Experiment Station, 
University of Missouri. Missouri Agricultural Experiment Station Journal Series 940. 

2 Stadler, L. J., Proc. Sixth Int. Cong. Genetics, 1, 274-294 (1932). 

3 Stadler, L. J., Cold Spring Harbor Symposia, 9, 168-177 (1941). 

4 Demerec, M., these PROCEEDINGS, 20, 28-31 (1934). 


5 Rhoades, M. M., Genetics, 23, 377-395 (1938). 
5 Rhoades, M. M., Cold Spring Harbor Symposia, 9, 138-144 (1941). 
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GRAPHICAL METHOD OF FACTORING THE CORRELATION 
MATRIX 


By L. L. THURSTONE 
THE PSYCHOMETRIC LABORATORY, THE UNIVERSITY OF CHICAGO 
Communicated April 14, 1944 


Multiple factor analysis starts with a square matrix R of order nm X n 
with cell entries ry, which are experimentally determined. The matrix R 
shows the correlation of each variable with every other variable in the test 
battery. The correlation matrix is symmetric since rx, = 7%,;. Some 
students write unity in the diagonal cells and this is legitimate for some 
problems. More often the diagonals contain the communalities h,? which 
are so chosen as to be consistent with the rank of the matrix as determined 
by the side entries. The factoring methods to be described can be applied 
for any set of diagonal entries.! 

The first object of multiple factor analysis is to find a factor matrix F of 
order m X r where r is the rank of the correlation matrix Ry. The factor 
matrix F must be written so as to satisfy the fundamental relation FF’ = R. 
The correlation matrix R can be factored into a factor matrix F in many 
different ways and F is not here unique. The current geometrical inter- 
pretation of R is that its entries 7, show the scalar products of pairs of test 
vectors J and K which are not necessarily of unit length unless this restric- 
tion is imposed by unit diagonals 7,;. Since the correlation coefficients are 
scalar products, it follows that no reference frame is implied by R. Each 
row j of the factor matrix F shows the entries a;,, which are the projections 
of the test vector J on a set of r orthogonal reference axes. The lack of 
uniqueness of F is geometrically interpreted in the free choice of an orthog- 
onal reference frame for the test configuration of ” test vectors J. The 
reference frame is explicitly defined in F relative to the test vectors but it is 
not defined in R. 

One factorial solution for F is to locate the first reference axis so as to 
maximize the sum of squares of test vector projections. The second refer- 
ence axis can be located by the same criterion applied to first factor correla- 
tion residuals, and so on until a set of r orthogonal axes have been deter- 
mined which correspond to rank r of the correlation matrix. In the writer’s 
original formulation of this problem, it was called the principal axes solu- 
tion.2 In the first paper on this problem a computational method was 
described for finding the principal axes solution F. It involved first factor- 
ing the correlation matrix into a factor matrix F with any arbitrary orthog- 
onal reference frame and a subsequent orthogonal transformation repre- 
senting a rigid rotation to the principal axes. That solution involved the 
characteristic equation of the matrix F’F and determination of its latent 
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roots. The next year Hotelling described an iterative method of factoring 
the correlation matrix directly into the principal axes solution which he re- 
named principal components.* Since the iterative methods, as well as all 
other available methods, of obtaining the principal axes solution are very 
laborious, it is of interest to students of multiple factor analysis to find 
shorter methods. 

Consider the first column of F for the principal axes solution. Its ele- 
ments are da;. It is desired to find the values of a; so that the first factor 
shall account for a maximum of the variance in the correlation matrix R. 
Then 


Vix — Ajdyi = Pixs (1) 


where aj; and a, are the first factor saturations of tests j and k, and p;, is 
the residual to be minimized. Squaring and summing for all correlation 
coefficients, we have 
D2? jx a 2220 py? jx LL j170y1” = Lr pix” == Z, (2) 
jk jk jk jk 
where the sum of the squares of the correlation residuals p;, may be denoted 
z. The partial derivative with respect to a; is then 


Oz ‘ss i 
en = Staten _ 4an2an (3) 
J 


from which we get 


An2ay* = IrpAn. (4) 

k k 
Let the sum Ya? = t. The summation Dra can be written as a matrix 
product Ra. Then 


Ra = at (5) 


where R is the given square matrix, a is a column vector, and ¢ is a scalar, 
namely, the largest latent root. The iterative process consists in starting 
with a trial vector u so that 


Ru = vru? (6) 


where (v=u*), or numbers proportional to these, constitute the trial vector u 
for the next iteration. The process continues until v2u? becomes propor- 
tional to the trial values u at which time the values of v (= u) are the de- 
sired values of a, and the first column of F is then determined. The num- 
ber of iterations can be considerably reduced by starting with R?, or R4, 
or even R®. While that looks effective in a text book example, it is not so 
inviting a procedure when we are confronted with a correlation matrix of 
order 70 X 70. 
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Consider the matrix multiplication Ru in equation (6) and, in particular, 
the scalar product of the row j of Ry, and the column vector u. This 
product is 


Ryu = v2? (7) 
and this can be written briefly as 


Ryu 


B42 
at 





vy = (8) 
but this is also the well-known equation for the slope of the regression line 
r on u through the origin for a plot of r against u. If we plot columns of 
rj, against the trial values u,, we get a plot with 7 points. The best fitting 
straight line through the origin (regression r on 7) can be drawn by inspec- 
tion and the slope of the line is easily read graphically with sufficient 
accuracy for the first few iterations. That slope is the value v;. To plot, 
say, 30 points and to draw by inspection the best fitting straight line 
through the origin takes less time than to add cumulatively 30 products on 
a calculating machine. 

By equations (5) and (6) it is evident that at the solution the values 
vj = Uj = a; When a set of values of v has been determined graphically, 
we plot v against « and find the slope m. In the first few graphical itera- 
tions, this slope will not be unity because the trial values u are probably too 
large or too small. If the true values are better represented by ku, in which 
k is a stretching factor, then a plot of r against ku will give a slope of v/k 
instead of v. When these graphs have been drawn, the slope of v/k against 
ku should be unity. But the obtained slope of v against u is m for the values 
actually used. Hence the observed slope m = k*?. Having found the slope 
m of the plot of v against “, we find k = ./m. Then if we should take a 
new set of trial numbers x; = ku, the iteration would give, instead of equa- 
tion (8), 


Rx 


Sy2 
ax 





= ¥, (9) 


where y = v/k. The slope of y against x should now be unity and the 
values of y should be used as the trial vector for the next graphical iteration. 
The determination of m and k can be done with a slide rule and the new 
trial values y = v/k can also be determined by the slide rule. The slope m 
can be found either by inspection or by simple summing with the method 
of averages in which m = 2v/Zw for like signed pairs of v and wu. 

The graphical method here described for the principal axes solution has 
also been adapted for a centroid method of factoring the correlation matrix 
by which the computations for a factor matrix can be done by a computer 
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without excessive labor even for a large matrix and without the use of 
tabulating machine equipment. That procedure is considerably faster 
and it will be described in a subsequent paper. 


A numerical example is here given. In table 1 we have a correlation 
matrix R of order 10 X 10 and rank 3.. The corresponding factor matrix F 
for the principal axes is shown in the same table. It is the solution to be 
found from R. The latent roots of R are 2.85, 1.47 and 0.63. These are 
the sums of squares of columns of F. In table 2 the initial trial vector u 
for the major principal factor is taken roughly proportional to the sums of 
columns of R. Column 1 of R is plotted against column u and the slope is 
estimated to be approximately +0.05 and it is recorded in column 2. 
Ten such plots give the values in column 7. The slope is read directly 
from each graph by noting the ordinate of a straight line fit at wu, = 1. 
Next obtain the two column sums Du, and 2v;. The ratio Sv/Su = m and 
k = Vm. This summation can be absolute sums for like signed pairs of u 
andv. Then compute y, = v,/k witha slide rule. These are also the trial 
values “2 for the next trial. Proceed likewise for three trials which give the 
desired values of a; to two decimals. 


The first factor residuals are then computed. These are py, = 7%. — 
GjQ,;1. Choose as a starting vector that column of the first factor residuals 
which has the largest absolute sum, ignoring diagonals. This is column 1. 
Its values are recorded in column x for factor two. The procedure is now 
the same as before. The third iteration was taken on a calculating machine 
and it gave the second factor loadings aj, the second column of F, to two 
decimals. The last iterations for each factor can be done on a calculating 
machine to obtain greater accuracy while the first few iterations can be done 
graphically to save time and labor. The computation for the third 
column of F was done in the same manner. 


Occasionally, when two or three latent roots are nearly the same, the 
iterative process will be found to oscillate or the convergence will be slow. 
This is an indication that the test configuration has nearly equal thick- 
nesses in these dimensions, so that the configuration is nearly circular or 
spherical in these dimensions. If the purpose is to extract the maximum 
variance from the correlation matrix, it does not matter where we place a 
set of orthogonal axes in those dimensions which are represented by nearly 
equal latent roots. In such a situation it is desirable to reduce the obtained 
vector v to a unit vector in the system without waiting for complete con- 
vergence. That can be done with a stretching factor p on v so that pu = a, 
where a are the desired factor loadings. The value of p can be found from 
the relation 1/p = V Suv. By this device the computer need not be 
unduly delayed by slow convergence in obtaining a factorially useful solu- 
tion. 
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The graphical method here described for estimating inner products can 
be applied in a variety of computational problems. 

1 The author wishes to express his indebtedness to the Carnegie Corporation of 
New York for a research grant in support of our development of multiple factor analysis 
and its application to the study of primary mental abilities. 

2 L. L. Thurstone, ‘Theory of Multiple Factors,”’ 1932, pp. 17-27, Edward Brothers, 
Ann Arbor, Mich. 

3 Harold Hotelling, “Analysis of a complex of statistical variables into principle com- 
ponents,” J. Educ. Psych., 24, 417-441, 498-520, (1933). 


THE DEVELOPMENT OF NORMAL AND HOMOZYGOUS BRACHY 
(T/T) MOUSE EMBRYOS IN THE EXTRAEMBRYONIC COELOM 
OF THE CHICK*t 


By S. GLUECKSOHN-SCHOENHEIMER* 
DEPARTMENT OF ZOOLOGY, COLUMBIA UNIVERSITY 


Communicated April 15, 1944 


The study of the causal morphology of mammalian embryos in spite 
of the great interest attached to it has not progressed very far because of 
the technical difficulties with which any experimental approach to the 
problem has met. Waddington and Waterman! grew rabbit embryos in 
vitro for a limited time and Nicholas and Rudnick* devised a method for 
raising rat embryos in a culture medium. In spite of these attempts, it 
was not possible to operate on mammalian embryos in early stages and 
let them continue their development inside the uterus, nor could such 
embryos be raised and develop normally in a suitable extra-uterine me- 
dium for any length of time. Methods which would accomplish this would 
be of importance not only for the study of normal causal morphology of 
mammalian embryos, but also for an experimental study of the embryog- 
eny of certain hereditary abnormalities. Searching for such a method, 
a new procedure for raising entire mouse embryos outside the uterus was 
developed and described in detail (Gluecksohn-Schoenheimer).* It con- 
sisted in removing the embryos from the uterus of the mother and trans- 
planting them into the extra-embryonic coelom of the chick embryo where 
they remained and developed for one or several days. This method has 
been applied in the experiments reported here to the study of embryos 
homozygous for the Brachyury mutation TT. 

As described by Chesley‘ the homozygous mutants (7/T) show extreme 
morphological abnormalities from the age of about 8 days on and die at 
about 10 days. The posterior body. region, including posterior limb buds, 
is missing completely and extensive abnormalities are found in notochord, 
neural tube and somites. 
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An additional abnormality, not mentioned by Chesley, was observed 
in the allantois. While the normal allantois, after having budded off from 
the embryo, appears turgid and proceeds to grow and to lengthen until it 
comes in contact with the ectoplacental cone, the abnormal allantois of the 
T/T homozygotes stays very short and has a shrunken tuberous ap- 
pearance. As a result of the abnormality of the allantois, umbilical ves- 
sels do not develop; thus the connection of the embryonic and the maternal 
circulation fails to be established and the embryo dies at the time when 
maternal circulation becomes all-important for the nutrition of the em- 
bryo (age of 10 days). 

Ephrussi,° in an attempt to determine the potentialities of the different 
tissues of the inviable 7/T embryos, explanted these tissues and grew them 
in tissue culture. Under the experimental conditions, the tissues sur- 
vived, grew and differentiated beyond the time and stage at which they 
would have died in utero. From his results, Ephrussi concludes that the 
death of the embryos in utero is probably caused by a “‘disturbance of 
the normal correlations’ and is due “exclusively to internal factors.”’ 
From our observations, it seems that the death of the embryo at the 
particular time of 10 days after fertilization is due to the failure of the 
circulatory connection with the mother rather than to any unknown internal 
factors. 

It seemed of considerable interest to study the extra-uterine development 
of whole 7/T embryos from two points of view: first, it was planned to 
examine their development in the chick host and compare it with that of 
embryos of the same genotype (7/T) in their normal environment; sec- 
ondly, it was hoped that by this approach we might be able to learn some- 
thing about the rdéle of the allantois in the development of the homozygotes. 

For this purpose, embryos derived from matings of heterozygous Brachy 
females (T/+) to heterozygous Brachy males (T/+) were used. A total 
of 149 transplants was made, 84 of which were successful. A transplant 
was considered successful if development of the mouse embryo proceeded 
after transplantation and if the embryo was alive when removed from the 
chick host. The mouse embryos were left in the chick embryo for 18 to 47 / 
hours. After this time, they were taken out from the chick host, exam- 
ined, recorded, fixed in Bouin, imbedded in paraffine, sectioned at 6 u 
and stained in Delafield’s Haematoxylin and eosin. The results to be 
reported here were arrived at after examination of the sections. 

As described before (Gluecksohn-Schoenheimer),* normal mouse em- 
bryos removed from the uterus of the mother 7 days after copulation (egg 
cylinder stage, Fig. 1) and transplanted into the extra-embryonic coelom 
of the chick develop and differentiate head structures, neural folds, 
somites, a beating heart and a tail bud (Fig. 2). The embryos from 
matings of Brachy (J7/+) by Brachy (T/+), transplanted into the extra- 











FIGURES 1-4 

Figure 1. Sagittal section through an egg cylinder stage of a normal mouse embryo, 
7 days 19 hours after copulation, = stage 11. A= Allantois. 80X. Figure 2. 
Sagittal section through a normal mouse embryo (No. 207) explanted into the extra- 
embryonic coelom of the chick in the egg cylinder stage (stage 11) and removed from 
there 49 hourslater. H = head,h = heart. 80X. Figure3. Sagittal section through 
a mouse embryo (No. 314) from a mating of T/+ X T/+, explanted into the extra- 
embryonic coelom of the chick in the egg cylinder stage (stage 11) and removed from 
there 24 hourslater. Embryo belongsinto ClassI = normal. H = head, S = somites. 
80X. Figure 4. Sagittal section through the same embryo (No. 314). Note turgid 
attached allantois. H = head, S = somites, A = allantois. 80x. 
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embryonic coelom of the chick, fall into three classes: the first class com- 
prises those embryos which develop normally in the host environment; 
they show good development of head, heart, somites and neural tube 
(Fig. 3); the allantois is big and turgid and attached to the ectoplacental 
cone (Fig. 4). The rest of the embryos show abnormal development and 
can, in turn, be separated into two classes, using the condition of the allan- 
tois as criterion. In both of these classes, II and III, the allantois is ab- 
normal, but in each class typically different from the other: Class II con- 
tains embryos with a short allantois which does not attach to the ecto- 





“A ALLANTCIS 


FIGURES 5-7 


Figure 5. Sagittal section through mouse embryo (No. 413) from mating of T/+ X 
T/+, grown in extra-embryonic coelom of the chick; Class II. No somites, under- 
development of head and trunk. Short allantois. 80%. Drawn from photograph. 
Figure 6. Sagittal section through mouse embryo (No. 312) from mating of T7/+ X 
T/+, grown ia extra-embryonic coelom of chick; Class III. No somites, underde- 
velopment of head and neural folds; note the short tuberous allantois. 80. Drawn 
from photograph. Figure 7. Photograph of a mouse embryo (7/T, Class III) which 
was explanted into the extra-embryonic coelom of the chick in the early head fold stage 
(8 days) and remained there for 47 hours. Typical features of T/T embryo: develop- 
ment of head, heart and anterior limb-buds, absence of posterior body structures. 
Compare Chesley‘, Fig.6. 20. 


placental cone; but while it is not big and turgid, it is not shrunken either. 
It looks as though it had been arrested in its development towards normal. 
The general condition and development of the embryos in this class are 
not as good as of those in the normal group: head and trunk structures are 
underdeveloped and somites are sometimes lacking entirely (Fig. 5). 
Whether there exists any causal relationship between the underdevelopment 
of the allantois and that of the rest of the embryo, or whether the abnormal 
allantois presents just one symptom of a general tendency to underdevelop- 
ment of the entire embryo cannot be decided now. The embryos of the 
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second abnormal group, Class III, are characterized by an allantois with 
a very typical abnormal appearance: it is small, not turgid but collapsed, 
and characteristically tuberous; it stays short and never attaches to the 
ectoplacental cone. Its most characteristic feature is its tuberosity (Fig. 
6). There is definite underdevelopment in the rest of the embryo; the head 
is small, the trunk short, the neural folds are thin, there are no somites, 
and hardly ever a beating heart. Four of the embryos in Class III de- 
veloped into typical T/T homozygotes (Fig. 7); all of these had been 
explanted in early head fold stages, i.e., at just 8 days after copulation. 
The rest of the embryos in this class which had been explanted in earlier 
stages—i.e., egg cylinder stage—could not be diagnosed by any of the 
typical features of 9-day T/T homozygotes because they did not develop 
far enough. But they were diagnosed by the typically abnormal condition 
of the allantois; and in addition they showed underdevelopment of head 
and trunk and absence of somites. 








TABLE 1 
EMBRYOS OBTAINED 
CLASS I CLASS II CLASS III 
NORMAL ABNORMAL ABNORMAL 7/T ? TOTAL 
Explantation stage of 
embryos, 
Stage: 7-10° 3 4 2 9 
Stage: 11-13” 8 22 11 5 46 
Stage: early head 
fold, 8 days 17 1 6 2 26 
Total 28 27 19 i 81 


* Early egg cylinder stages, end of 6 days. 
> Egg cylinder stages, 7 days, beginning growth of allantois, slight thickening of head 


folds. 


When we consider the numerical distribution of the explanted embryos 
among the three different classes (table 1), we find that 19 out of a total 
of 81 fall into Class III, i.e., the group comprising the embryos supposed to 
represent the abnormal 77 class. Of these 19, only 4 developed all of the 
typical TT features; but all of the others had the typically abnormal 
tuberous allantois which is found in all 77 embryos developing in their 
normal environment (see above) and which never appears in explanted 
embryos from matings of normal by normal. We expect one quarter of the 
embryos from T/+ by T/+ matings to be T/T and the actual number of 
19 out of 81 fits expectation very well. 

The second class of abnormals—27 embryos—might possibly represent 
the T7/+ embryos. If one considers the different explantation stages 
separately in this class, it appears that of the embryos explanted at earlier 
stages (7-13) about 50% fall into. Class II while of those explanted at 
later stages (early head fold and 8-days stage) about 70% fall into Class I 
and only one embryo is classified as Class II. It is conceivable that the 
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developmental potencies of the very early 7/+ embryos under the un- 
favorable explantation conditions are reduced in comparison with those of 
+/+ embryos. If this is the case, it might well be that the allantois, which 
is completely abnormal in-the homozygotes, reacts to the abnormal ex- 
plantation conditions more strongly in the young heterozygous embryo 
than it does in the homozygous normal, showing lower developmental po- 
tencies under abnormal conditions than does the normal. In this con- 
nection it seems of interest that a deficiency of the number of heterozygous 
short-tailed mice (7/+) appeared in matings of Brachy short-tailed 
mice (7'/+) inter se (Kobozieff and Pomriaskinsky-Kobozieff)* and in our 
own laboratory (Dunn).’ Although some of this deficiency could be 
accounted for by the failure of the Brachy mutation to express itself in 
some of the heterozygotes, a differential prenatal mortality of 7/+ 
embryos might be responsible for the lower number of heterozygotes in 
addition to the normal appearance of some 7/+ mice. The total num- 
ber of embryos in Class II does not correspond to the expected 50% of 
T/+ embryos while the number of 28 in the normal group is higher than 
the expected 25%. One would thus have to assume that a certain fraction 
of 7’/+ embryos, namely, the older ones, were able to develop normally, 
or one might ascribe the distribution of embryos among the two classes 
entirely to chance. In the latter case, one would assume identical develop- 
mental potencies for +/+ and T/+ embryos, and they would be equally 
distributed between Classes I and II. Their total of 55 corresponds well 
to the expected 75% of +/+ and T/+ embryos combined. While, thus, 
there remain some doubts about ascribing the embryos in Classes I and II 
to two genetically different groups, +/+ and 7/+, the remaining 25% 
fall into a group distinctly and characteristically different from the rest 
of the embryos, independent of the stage of explantation; this group is 
assumed to represent the 7/T class. 

It thus appears that the abnormal development of 7/7 embryos is 
determined before it is apparent morphologically and is independent of the 
stage of explantation (table 1). The earliest stages which it is possible to 
expiant, i.e., early egg cylinder stages, show the tendencies to abnormal 
development of the entire embryo. Earlier stages cannot be tested, since 
they do not develop at all under the explantation conditions described, 
regardless of whether they come from crosses of normals or Brachys. 

An analysis of the rdéle of the allantois in the development of the normal 
and abnormal embryos was not possible on the basis of our experiments. 
The fact that an abnormal or underdeveloped allantois occurs always in 
combination with abnormal or subnormal development of the entire 
embryo seems highly suggestive. The allantois of the mouse might well 
serve as an organ for intake of nutrient substances and excretion of waste 
products. Any disturbance in the function of this organ might thus lead 
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to secondary abnormalities throughout the entire embryo. However, 
experiments involving the extirpation of the allantois and studying its 
effects on the embryo would be required to settle the question of the réle 
of the allantois in the development of the normal mouse embryo. Such 
experiments might also provide an answer to the question of the possible 
causal réle of the allantois in the abnormal development of 77 embryos. 

Summary.—Embryos from crosses of Brachy short-tailed mice (7/+) 
inter se were removed from the mother in early stages of development and 
transplanted to the extra-embryonic coelom of the chick. They re- 
mained in the host for 18 to 47 hours and were then removed for study. 
The embryos obtained could be separated into three different classes: 
(1) normal embryos; (2) abnormal embryos with different degrees of ab- 
normalities, with reduced developmental potencies, possibly representing 
the 7/+ embryos; (3) abnormal embryos with typical abnormalities, 
representing 7/T embryos. 7/T embryos will thus develop true to their 
genotype, even if explanted before the time at which morphological differ- 
ences appear. The possible causal réle of the allantois in the development 
of these embryos is discussed. 

+ This research was aided by grants from the Josiah Macy, Jr., Foundation and the 
Fund for Research of Columbia University. 

* The author wishes to express her gratitude to Dr. L. C. Dunn for his suggestions 
and his constant interest during the course of this work. 
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REMARK ON A DISTORTION TENSOR FOR ELASTIC 
DISPLACEMENTS 


By T. Y. THomas 


UNIVERSITY OF CALIFORNIA AT LOS ANGELES AND INSTITUTE FOR ADVANCED STUDY 


Communicated April 18, 1944 


One of the most successful of the various conditions for yield resulting 
from a state of stress produced by mechanical means is that due to von 
Mises.' This condition states that yield will occur when the energy per 
unit volume produced by a change of shape alone reaches a certain value 
characteristic of the material under consideration. It would thus appear 
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that distortion, i.e., change of shape as distinct from change of volume, is 
one of the primary factors in the occurrence of yield. This suggests that in 
addition to the ordinary tensors of stress and strain (related by Hooke’s 
law) we consider a distortion tensor and seek to express the condition of 
yield directly in terms of this tensor. 

A differentiable point transformation x — Z referred to rectangular co- 
ordinates and having its jacobian everywhere different from zero (and 
hence locally non-singular) will be conformal (i.e., angles will be left un- 
altered) if, and only if, 


af Oz* OF" 
2 Pony — Phen 


where p is some positive function of position. When the transformation is 
an elastic displacement it is customarily represented by @* = x* + u*(x) in 
which the u* (also denoted by u, since rectangular coérdinates are em- 
ployed) are the components of a vector. In accordance with the usual 
assumptions, the u’s and their derivatives are small quantities, the product 
of which can be neglected. Hence 

= = 06 + u%,, (2) 


ox? 


(p > 0), (1) 


and 
5s, + Ug,y + up = p5p,, (3) 


when (2) is substituted into (1). Taking the determinant of both members 
of (1) and making use of the form (2) of the derivatives 07*/0x* we have 


\06 + u%q|? = |pdag| = 0°, 
or 
p= |ig tum = 1+ %/u% =1+%sd, (A= u%). 
Substitution of this value of p into (3) gives 


cup ~ = 8p = — (Cap = "Valtass + Mp, 0)) (4) 


where e is the strain tensor of elasticity theory. The quantity A is the 
dilatation or increase per unit volume due to the displacement. Conversely 
if (4) is satisfied the elastic displacement is conformal. 

We now define a tensor D as an invariant of an arbitrary elastic displace- 
ment by putting 


A 
Dap = Cap — 3 Sap: 
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We call D the distortion tensor. From the above, the displacement will be 
conformal (i.e., without distortion) if, and only if, D = 0. We can there- 
fore think of the tensor D, in general, as giving a measure of the distortion 
produced by the elastic displacement.’ 

It is well known that any elastic displacement causes a displacement of 
the local neighborhood of an arbitrary point P which is the result of the 
following three displacements, namely (a) a translation P — Q, (8) a rigid 
rotation about Q, and (y) a non-rigid displacement or deformation having Q 
as a fixed point. This latter deformation can be resolved, on the basis of 
the distortion tensor D, into a component with distortion and a component 
without distortion. Thus if 6x denotes an arbitrary small vector issuing 
from the above point Q, the part of the deformation due to the distortion 
will be represented by D,,5x® and the part without distortion by 
(A/3)éqg5x* these quantities being the components of the changes pro- 
duced in 6x* in this decomposition of the deformation. 

Hooke’s law for a homogeneous and isotropic medium can now be written 


Eup = 2uD 28; 6 = (3d > 2) A, (5) 
where 
| a = Tap =" 5 bap 8 = br oe 


Thus E£ is the distortion analogue of the stress tensor o while \ and uv are the 
usual elastic constants. Now denote by 7, and y, (7 = 1, 2, 3) the principal 
values (at any point) of the tensors o and E, respectively. From the 
definition of these quantities we immediately have the relations y; = 7,;— 
6/3 and when written in full these become 


vm = Y/s[(71 — 72) + (11 — 78)], | 
12 = 1/3[(T2 iia T3) a (m1 7 T2) |; ( (6) 


vs = “/al—(n1 — 2) — (m2 — 19)]: 
Now consider the square of the tensor E. Denoting this by E* we have 


E* = LE .pE op = 2 (Yabas)(Ya5ae) = 21% 
1/s[(71 ore T2)? + (71 sa 73)” + (rn — T3)*], 


when use is made of the above relations (6). But apart from a constant 
factor the above bracket expression is the energy per unit volume due to 
distortion.* Since the squares of the tensors D and E differ only by a con- 
stant factor on account of (5), the von Mises yield condition can therefore 
be stated as follows: yield will occur when the square of the distortion tensor 
D (or its stress analogue E) attains a characteristic value. 

By a distortion invariant of an elastic displacement we understand any 
invariant which depends solely on the distortion. In particular any scalar 
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function of the components ogg (or é,g) will be a distortion invariant if, 
and only if, it is expressible as a function of the components E, (or D,,). 
The sufficiency of this condition is immediately obvious. To show the 
necessity for a scalar S(o,g) we have 


S(Gap) = S(tadap) = Sl(ta + 7)dag] 


where 7 is arbitrary. But choosing r = —0/3 the latter function is seen to 
become S(y.5.8)0r S(E,s) from which the conclusion follows. This 
result can obviously be extended to differential and non-scalar invariants. 
The above quantity E* or D? is an example of a distortion scalar. So like- 
wise is the quantity y1 — ys = 71 — 73 which is proportional to the maxi- 
mum shearing stress (when 7; 2 72 2 73). This stress also appears to be of 
significance in determining the condition of yield.‘ Other examples are 
the magnitudes of the principal values y; of the tensor E (which differ from 
the principal values of the distortion tensor by a constant factor). Useful 
modifications of the known yield conditions may possibly be given by these 
latter or other distortion scalars which can be constructed. Such invari- 
ants will automatically take care of the failure to reach a state of yield as 
the result of hydrostatic pressure or uniform tension. 


1 Nachr. d. Gesellsch. d. Wissensch. zu Géttingen, Math.-Phys. Klasse, pp. 582-592 
(1913). See also, Nadai, A., Plasticity, McGraw-Hill, Chap. 13 (1931), and Timo- 
shenko, S., Theory of Elasticity, McGraw-Hill (1934), p. 137. 

2 The usual formula for the change 50 of an angle 6 due to an elastic displacement gives 
this change in terms of the components égg of the strain tensor. Cf. Love, The Mathe- 
matical Theory of Elasticity, Cambridge, Vol. 4 (1927), p. 62. Also, Timoshenko, §&., 
loc. cit., p. 192. By making use of the above definition of the components Dag it can 
be shown that 60 is given by 


sin 0 80 = (D; + Dz) cos 0 — 2Dapr%4 
where, 
D, = Dapvth, Da = Dapo, 
and »; and v2 are unit vectors determining the angle 6. Hence the change in an angle 
produced by an elastic displacement is expressible directly in terms of the distortion 
tensor D. 


3 See Nadai, S., loc. cit., p. 72, and Timoshenko, &., loc. cit., p. 137. 
4 See, ‘or example, Nadai, A., loc. cit., p. 60. 
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THE MINIMUM RANK OF A CORRELATION MATRIX 
By A. A. ALBERT 
DEPARTMENT OF MATHEMATICS, UNIVERSITY OF CHICAGO 
Communicated April 22, 1944 


A correlation matrix is a symmetric matrix R with prescribed non- 
diagonal elements (correlation coefficients) and the matrix problem with 
which we are concerned is that of determining the diagonal elements (com- 
munalities) of R so that the resulting R will have minimum rank p. Ina 
recent article! I defined the ideal rank r of R to be the largest order of a non- 
vanishing minor of R such that all elements of the minor are correlation 
coefficients. Evidently p 2 r and I gave a procedure for determining the 
communalities under the hypothesis that p = r. 

Indeed we may always select our notation for the elements of M so that 


A B'G’ 
R={|{BC H'}, 
GHK 


where A, B, C are r-rowed square matrices and the matrix B has an inverse 
B-'=E. Then I proved the 

Lemma 1: [If the ideal rank of Ris r the matrix K — GEH’ is a diagonal 
matrix. The rank of R is r if and only if K = GEH', H = GEC, G = 
HE’A, B’ = AEC. 

In cases usually encountered in factor analysis the ideal rank r is much 
less than the order » of R and the r-columned matrices G and H normally 
have their maximum rank r. The purpose of the present note is to provide 
a proof of the rather astonishing 

THEOREM. Let the ideal rank r of R be also the rank of its submatrices G 
and H. Then there exist unique communalities such that the resulting correla- 
tion matrix R has rank r. 

We shall prove this result by several applications of the first statement in 
Lemma 1. Since K is symmetric it implies that GEH’ is symmetric. Let 
the rank of G be r so that we may select a notation for R so that 


6) 1-0) 


for r-rowed square matrices S and T such that Sis non-singular. Define a 
non-singular matrix 


eS I ’) 
~ \—PS-1I 
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of n — 2r rows, for appropriately selected identity matrices J, and form the 


products 
fs eo 
LG = a LH = al 


where W = Q — S-! PT. Then 


L(GEH’)L' = (LG)E(LH)’ = oy (T’W’) = ors ge 


is a symmetric matrix and necessarily SEW’ = 0. Since SE is non- 
singular we have proved that W’ = 0, W = 0, 


Me i enc: piel Se\., . (ars 
LH = (5) L(GEH’)L' = L post PAB = ( : = 


Suppose now that H has rank r so that LH has rank r and T must have 
rankr. Then T is non-singular and we may prove 

LemMa 2: The matrices A — B'T—!S and C — BS~'T are diagonal 
matrices. 

For the submatrix of R consisting of its first 3r rows and columns is a 


correlation matrix 
A B'S’ 
N =O PE 4 
ae ley Bi 


which is ::xactly like R except that K has been replaced by U, G by S, H by 
T. Every diagonal element of N is a diagonal element of M and B is an r- 
rowed non-zero minor of N. Then the ideal rank of N is r and, by the first 
part of Lemma 1, U — SB~'T”’ is a diagonal matrix. However, this result 
is actually a part of the statement that K — GEH’ is diagonal. We now 
obtain two correlation matrices 


UTS A S’ B’ 
M=({T'C BL M=|(SUT 
S’ B’ A BT’C 


of ideal rank r from N. The matrix JN, is obtained by interchanging the 
first and third blocks of r rows and columns in N and the matrix N2 by 
interchanging the second and third blocks of r rows and columns in N. 
Our lemma then follows from the first part of Lemma 1. 

If R has rank r then N, Ni, Nohaverankr. Applying Lemma 1 we have 
the uniqueness part of our theorem which we state as 

LemMa 3: Let Rhaverankr. Then the communalities are uniquely deter- 
mined by the equations 
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A = B'T"S, C = BST, K = GBH’. (1) 


This lemma provides a modification of my earlier technique for comput- 
ing the communalities. In the new technique we compute the inverses of 
the three matrices S, T, B and the diagonal elements of three matrix prod- 
ucts. 

To complete the proof of our theorem we assume now that R has ideal 
rank r and that 7 is the rank of Gand H. By Lemmas 1 and 2 we may 
select the diagonal elements of R so that (1) is satisfied. Then 


L(H — GEC) = oe - +c, = ia nie =0 


since C = BS-'T, T = SEC. Similarly 


L(G — HE’A) = (3) = Po oe (° - pas — 


since A = B’/T—'S, S = TE’A. But L is non-singular and therefore H = 
GEC, G = HE’A. Finally, we have the identity 


GE(B’ — AEC)’ = (H — GEC)E’A + G — HE’A 


so that, since the right member is zero, GE(B — AEC)’ = 0. The first 
r rows of this product form the matrix SE(B’ — AEC)’ = 0. Thus 
(B’ — AEC)’ = 0, B’ = AEC, the rank of Ris r. This completes the 
proof of our theorem. 

Let us close by observing that R cannot have its ideal rank as its rank 
unless the ranks of G and of H are the same. For this is a consequence of 
the fact that when H = GEC the rank of H is at most the rank of G, and 
that when G = HE’A the rank of G is at most the rank of H. However, 
even when G and H have equal ranks /ess than r the minimum rank of R 
need not be its idealrankr. Indeed let G = H = 0, K bea diagonal matrix 


and 
_ {10 _ (x 0 _ {nl 
B=(oi) 4-(c2) ¢= (75) 


Then every minor of R with more than two rows and containing no com- 
munality has a zero row and sor = 2. But B’ = AEC now becomes C = 
A~-‘ and this is false since A— is a diagonal matrix and C is not a diagonal 
matrix. 


1 The matrices of factor analysis, Proc. Nat. Acad. Sci., 30, 90-95 (1944). 





